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Abstract 



We apply the realistic shell model which includes the coupling between many-particle 
(quasi-)bound states and the continuum of one-particle scattering states to the spec- 
troscopy of mirror nuclei: 8 B and 8 Li, as well as to the description of low energy 
cross sections (the astrophysical 5 factors) in the capture reactions: 7 Be(p, 7) 8 B and 
7 Li(n,7) 8 Li . 
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I. INTRODUCTION 



The solution of solar neutrino problem, i.e., an observed deficit of neutrinos with respect 
to predictions of the Standard Solar Model (SSM) [[| , is passing through an understanding 
of the capture reaction: 7 Be(p, 7) 8 B. The 8 B produced in the solar interior in the reaction 
7 Be(p, 7) 8 B, is the principal source of high energy neutrinos detected in solar neutrino ex- 
periments. Therefore, if we stay within the framework of standard electroweak theory and 
the SSM 0, the observed deficit of those neutrinos can be traced back, at least partially, 
to the value of low energy 7 Be(p, 7) 8 B radiative capture cross section which determines the 
magnitude of neutrino flux and remains the most uncertain quantity in the SSM 0. At the 
solar energies ( ~ 20keV), this cross-section is too small to be directly measurable. For this 
reason, the theoretical analysis of this reaction is so important. On the other hand, whenever 
the measurement is feasible {Ecu > 150 keV), the exact value of the capture cross section 
depends : (i) on the normalization obtained indirectly from the 7 Li(<i,p) 8 Li cross section 
and, (ii) on the model dependent extrapolation of measured values of the cross-section down 
to the interesting domain of solar energies. Experimental data for the 7 Be(p, 7) 8 B capture 
cross section are varying strongly [0-0] , though more recent experiments consistently indi- 
cate low values (S < 20 eV-b) for the astrophysical factor S = <jcm{Ecm)Ecm exp(— 2^77), 
where r\ = e 2 ZiZ 2 /hv [§|,[7j. Also the Coulomb dissociation experiments deduce low value 
for the astrophysical S - factor though this value depends strongly on the amount of E2 - 
contribution in the Coulomb dissociation which is not yet completely understood. 

Many different theoretical models have been applied for the calculation of S - fac- 
tor at stellar energies PHl6ll, and their predictions are often in striking disagreement 
among each other, confirming strong model and/or approximation dependence of calcu- 
lated cross-section. Part of the theoretical ambiguities can be removed by a simultaneous 
study of the 7 Li(n, 7) 8 Li mirror reaction, which has also been studied by several experi- 
mental groups |T7|-|l9|]. Also for this reaction, largely different values for the low energy 
capture cross section have been reported [117- T^]. In the context of the solar neutrino 
problem, the 7 Li(n, 7) 8 Li cross section is often used to extrapolate the capture cross sec- 
tion for the reaction 7 Be(p, 7) 8 B down to the solar energies at E C m ~ 20keV || . But 
the 7 Li(n, 7) 8 Li reaction cross section at very low energies is also extremely interesting 
in itself as it provides the essential element of rapid process of primordial nucleosynthe- 
sis of nuclei with A > 12 in the inhomogeneous big-bang models [f2~0|-f23|| . Indeed, in 
the inhomogeneous big-bang hypothesis [p0|-|23|, the main reaction chain leading to the 
synthesis of heavy elements is [^] 1 H(n, 7) 2 H(n, 7) 3 H(<i, n) 4 He(t, 7) 7 Li(n, 7) 8 Li, and then 
8 Li(a, n) 11 B(ra, 7) 12 B( / 9~) 12 C(n, 7) 13 C, etc., for heavier nuclei. In this sense, the reaction 
7 Li(n, 7) 8 Li is a key process to bridge the gap of mass A = 8 and to produce heavy ele- 
ments. 

The theoretical description of weakly bound exotic nuclei close to the drip-line, such as, 
e.g., 8 B on the proton rich side of the drip line or 11 Be, n Li on the neutron rich side of 
the drip line, is challenging due to the proximity of particle continuum which implies the 
strong modification of effective nucleon - nucleon interaction and causes the unusual spatial 
properties (halo structures, a large diffusivity) of nucleon density distribution . Nowadays, 
these properties are perhaps somewhat better understood near the neutron drip line than 
near the proton drip line. In weakly bound exotic nuclei, number of excited bound states or 
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narrow resonances is small and , moreover, they couple strongly to the particle continuum. 
Hence, these systems should be described in the quantum open system formalism which does 
not artificially separate the subspaces of (quasi-) bound (the Q-subspace) and scattering 
(the P-subspace) states. For well bound nuclei close to the (3 - stability line, microscopic 
description of states in the first subspace is given by the nuclear shell model (SM) with model- 
space dependent effective two-body interactions, whereas the latter subspace is treated in 
terms of the coupled channels equations. For weakly bound exotic nuclei, the validity of 
this basic paradigm is certainly questionable, and we propose to change it by considering 
the approximation which takes into account coupling between Q and P subspaces in terms 
of the residual nucleon-nucleon interaction. This coupling will consistently modify both the 
scattering solutions and the spectroscopic quantities for interior bound states. 

One possibility for such an approach could be the Continuum Shell Model (CSM), which 
in the restricted space of configurations generated using the finite-depth potential has been 
studied for the giant resonances and for the radiative capture reactions probing the micro- 
scopic structure of these resonances p4|-f27|j . This approach may be insufficient for the 
description of low lying excitations, in particular in the nuclei close to drip lines, where it is 
essential to have a most realistic description of bound states. For that reason, the corner- 
stone of our approach, which is called in the following the Shell Model Embedded in the 
Continuum (SMEC), is the realistic SM itself which is used to generate the iV-particle wave 
functions. This choice implies that coupling between the SM states and the one-particle 
scattering continuum must be given by the residual nucleon - nucleon interaction. The first 



application of the SMEC model has been published recently in Ref. [29[ . 

As said before, we are interested in describing low lying bound and resonance states in 
exotic nuclei. For that reason, we restrict the description of particle continuum to the subset 
of one-nucleon decay channels. This should be a reasonable starting point for describing 
both the microscopic structure of 8 B and 8 Li and the corresponding reactions: 7 Be(p, 7) 8 B 
and 7 Li(ra, 7) 8 Li. At higher energies, e.g., above a or t emission thresholds, the one-particle 
continuum approximation prohibits a reliable description of more complicated multi-nucleon 
decay channels as well as the residual correlations generated in the many-body wave functions 
of bound states by the coupling to those channels. In principle, this obstacle can be removed 
in future studies and further improvement of the SMEC (CSM) to include more complicated 
channels like, e.g., the a - decay channel, can be done following the approach of Balashov et 



al. |30[ . This rather involved extension of the SMEC will not be discussed here any further. 
One should also be aware that in the case of two-nucleon halo nuclei or at higher energies 
above the threshold for the three-particle decay, the one-particle continuum approximation 
is oversimplified. 

The paper is constructed in the following way. In Sect. II, we present some details of 
the SMEC model, stressing in particular those elements which differ this model from the 
CSM. Sect. Ill is devoted to the discussion of spectroscopic properties of 8 B (Sect. III.B) 
and 8 Li (Sect. III.C). We shall discuss influence of the residual interaction which couples 
Q and P subspaces on the energy and the width of many-body states. We will also discuss 
an important self-consistency correction to the average single particle (s.p.) potential which 
results from this residual interaction. The self-consistent average potential is then used to 
generate radial formfactors of many-body states which enter the coupling matrix elements 
between Q and P subspaces. Sect. IV is devoted to the discussion of 7 Be(p, 7) 8 B and 
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7 Li(n, 7) 8 Li capture cross sections. Finally, the most important conclusions are summarized 
in Sect. V. 



II. THE SHELL MODEL EMBEDDED IN THE CONTINUUM (SMEC) 

In the decaying nucleus or in the reaction processes such as, e.g., the radiative capture 
process, all asymptotic channels are given by the solutions with outgoing waves. Such 
physical systems, in which bound (localized) interior states are coupled to the asymptotic 
scattering channels, are called the open systems. One way out of this complication is to 
describe the quantum open system together with its environment of asymptotic channels as 
the quantum closed system. This implies that the projection operator technique is used to 
separate P subspace of asymptotic channels from the Q subspace of many-body states which 
are build up by the bound s.p. wave functions and by the s.p. resonance wave functions. 
In this latter case, one should use certain cutoff radius to define this part of resonance 
wave function which belongs to the Q subspace and the remaining part is put into the P 
subspace. In our case, the P subspace contains (N — 1) - particle states with nucleons on 
bound orbits and one nucleon in the scattering state, but also this part of the s.p. resonance 
wave functions which is outside of the cutoff radius. In fact, we define the P subspace from 
the condition 

P + Q = i , (i) 

which also implies that all wave functions used in the calculations are orthonormalized in 
the usual sense. Physically, the closeness assumption for the total system means that the 
states of (N — 1) - nucleus which define the asymptotic channels for the studied reaction are 
assumed to be stable. This assumption in turn means that the most important channels are 
supposed to be those which involve low-lying states of the residual nucleus. In the domain 
of low energy excitations and/or for weakly bound nuclei away from the /3-stability line this 
is most likely a good approximation. However, in the domain of higher excitation like, e.g., 
for the giant resonances, the quality of this approximation which depends on the width of 
considered many-body states in N — 1 nucleus, may be hazardous for broad states. 

The key element of the CSM is the treatment of s.p. resonances, which on one side may 
have an important amplitude inside a nucleus and, on other side, they exhibit asymptotic 



behaviour of scattering wave functions Pj| . The part of resonance for r < R cut , where 
R cut is the cut-off radius, is included in Q subspace, whereas the remaining part is left in 
the P subspace pi . The wave functions of both subspaces are then renormalized in order 



to ensure the orthogonality of wave functions in both subspaces. It should be mentioned 



that in the earlier formulation of the CSM | 28fl , the contribution of the s.p. resonances both 
to Q and P subspaces have been neglected. 



In the SMEC calculations, we solve similar equations as in the CSM |24j but , as explained 
below, due to specificity of exotic nuclei, we change certain important ingredients of the 
CSM. For the bound states we solve the standard SM problem: 

H QQ $ i = E i <S> i . (2) 

Hqq = QHQ is identified with the SM Hamiltonian and its $j are the iV-particle (quasi-) 
bound wave functions. We believe, that for a quantitative description of low lying states in 
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the exotic nuclei one has to use as a starting point the accurate many-body wave functions 
in the Q subspace which are provided by the SM with effective interactions. 
The residual coupling of P and Q subspaces is given by the zero-range force 

Vl2 = -V$[a + (1 - a)Pf#( ri - r 2 ) , (3) 

where P° 2 is the spin exchange operator. We assume that the SM Hamiltonian Hqq contains 
already effects generated by this residual coupling and we do not modify Hqq anymore. 
For the continuum part, we solve the coupled channel equations: 

- #PP)d +) = , (4) 

where index c denotes different channels and Hpp = PHP . The superscript (+) in (|J) 
means that boundary conditions for outgoing scattering waves are used. In our case, we have 
ingoing wave in the input channel and outgoing waves in all other channels. The structure 
of (N — 1) - nucleus is given by the SM, whereas one nucleon occupies a scattering state. 
The channel states are defined by coupling one nucleon in the continuum to a many-body 
state of (N — 1) - nucleus. 

The SM wave function has an incorrect asymptotic radial behaviour for unbound states. 
Therefore, to generate both the radial s.p. wave functions in the Q space and the scattering 
wave functions in P space we use the finite-depth average potential of Saxon- Woods (SW) 
type with the spin-orbit part included: 

U(r) = Vof{r) + Vso(41-8)±yQ + Va , (5) 
where f(r) is the spherical symmetric SW formfactor: 

(6) 

The Coulomb potential Vq is calculated for the uniformly charged sphere with radius Rq. 
Coupled channel eqs. (Q) can be written more explicitly as: 

E(£ (+) - H c >)d {+) = , (7) 

c 

where 

H c > = (T + U)5 CC , + v J c , . (8) 

In the above equation, T stands for the kinetic-energy operator and v J , is the channel- 
channel coupling generated by the residual interaction. The explicit formula for v J cc , is given 
in flX4]). The s.p. potential in (||) consists of an 'initial guess' U(r), and the diagonal part of 
the coupling potential v J cc which depends on both the s.p. orbit (pij and the considered many- 
body state J n . Obviously, this correction cannot be neglected when generating s.p. wave 
function (pij for a given J % . These s.p. wave functions define Q subspace and thus modify 
the diagonal part of the residual force. So this implies a self-consistent iterative procedure, 
because the change of s.p. wave function changes in turn the correction coming from the 



f(r) 
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residual force (see the Appendix A). As long as this correction is small in comparison with 
the initial average potential, the iterative procedure is fastly converging to the new self- 
consistent average potential: 

l/(«0( r ) = Z7(r)+u^(r) . (9) 

The parameters of the initial average potential (|5|) are chosen in such a way that the re- 
sulting potential U^ sc \r) reproduces energies of experimental s.p. states, whenever their 
identification is possible. In the next section, we will show examples of such self-consistent 
potentials. The dependence of the generated correction in the potential @ on isospin of 
s.p. states generates a mutual dependence of average potentials for protons and neutrons in 
the iterative procedure. Using self-consistently determined radial wave functions instead of 
those generated by U(r) means that the matrix elements V a ^s (Al) of the residual force V\i 
depend not only on the s.p. wave functions <pij involved but also on the many-body state J 71 . 
This is another interesting aspect of the self-consistent procedure determining the average 
potential in SMEC. 

The third system of equations in CSM consists of inhomogeneous coupled channel equa- 
tions: 

(£(+) _ H PP )4 +) = HpQ$i = Wi (10) 

with the source term which is primarily given by the SM structure of N - particle wave 
function for state $j . The explicit form of the source term is given in (AT). These equations 
define functions uo\ +S> , which describe the decay of quasi-bound state $j in the continuum. 
The source Wi couples the wave function of N - nucleon localized states with (N — 1) - 
nucleon localized states plus one nucleon in the continuum. For zero-range residual force (0), 
formfactor of the source term is given by the s.p. wave functions of the same self-consistently 
determined average potential U^ sc '(r) as used to define the subspaces Q and P. 

The full solution of the many-body problem can be expressed by means of three distinct 
functions: , £jg and Ui , and reads : 



= C E + E( $ » + "0 „ TT eff < *J I H QP I C E > , (11) 



'■J 



E-m 



where 



Hqq = Hqq + HqpG^Hpq , (12) 



is the effective SM Hamiltonian including the coupling to the continuum, and is 
the Green function for the motion of s.p. in P subspace. Matrix Hqq is non-Hermitian 
(symmetric and complex matrix) for energies above the particle emission threshold and 
Hermitian (real) for lower energies. Consequently, the eigenvalues of Hqq are real for bound 

states and complex for decaying states. Matrix Hqq can be diagonalized by the orthogonal 
transformation: 
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with complex eigenvalues: Ei — , which depend on the energy E of particle in the 
continuum. The eigenvalues of Hqq at energies Ei(E) = E, determine the energies and 
widths of resonance states. With these changes, one obtains: 

V C E = £ C E + Y£li ^_ ^<$ i \H\&> (14) 

i E - Ei + (i/2)r< ' 1 £ v ; 

for the continuum many-body wave function projected on channel c , where 

roc | 

^-^ + E/ £c dE ^E' Ei+) _ E ' <& \n\*i> , (is) 

is the wave function of discrete state modified by the coupling to the continuum states. 
Hence, this formalism if fully symmetric in treating the continuum and bound state parts 



of the solution : ^> C E (in Eq. ([14]) ) represents the continuum state modified by the discrete 
states, and Oj (in Eq. flT5])) represents the discrete state modified by the coupling to the 
continuum. More informations about those features of SMEC which are the same as in the 
CSM, can be found in the paper by Bartz et al. [El]- 



III. SPECTROSCOPY OF 8 B AND 8 LI NUCLEI 

A. The self-consistent average potential spanning the Q subspace 

The essential element of SMEC approach is the construction of Q - subspace. As ex- 
plained in the previous section, this is obtained by an iterative procedure which for a given 
initial average s.p. potential (|) and for a given residual two-body interaction (|3|) yields the 
self-consistent s.p. potential depending on the s.p. wave function 0y, the total spin J of the 
iV-nucleon system as well as on the one-body matrix elements of (N — 1) - nucleon daughter 
system. The parameters of different initial SW potentials for [p ® 7 Be] and [n (g) 7 Li] systems, 
which are used in this work are summarized in Table |. All these potentials have the same 
parameters of radius Rq = 2.4 fm, surface diffuseness a = 0.52 fm, and spin-orbit coupling 



Vso — — 4MeV. Cohen - Kurath (CK) interaction [31] is used as a SM interaction. The 
potentials of Table | have been obtained for different parameters (1 — a) of the spin-exchange 
component in the residual interaction (|3|) and for different s.p. orbits <pij which correspond 
to energies £y (see the second column in Table |) in the self-consistently determined poten- 
tial U( sc \r). In all cases, the strength of the residual interaction (|3|) is V$ = 650MeV-fm 3 . 
This value yields a reasonable compromise between the description of 8 B and 8 Li spectra 
and in particular their decay widths. The dependence of the energy spectrum of 8 B on V$ 
will be presented in the Sect. III.B. 

In Fig. 1 , we show examples of calculated potentials in 8 B for the proton s.p. orbital 
l|?3/2 in two different total spin states: J 71 " = 2 + and J 71 " = 1 + . The calculations have been 
performed using the initial potential U(r) for (1 — a) = 0.27, which is chosen in such a 
way that the self-consistent potential U^ sc \r) yields proton s.p. orbit bound at the 

experimental binding energy of ground state (g.s.) J n = 2f. The same U(r) is then taken 
both for 2 + and 1 + states. The spectroscopic factor of proton lpz/2 s.p. state in the g.s., is 
close to 1 for CK interaction This allows to identify position of proton lpz/2 s.p. orbit 
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in J n = 2 + state, i.e., we demand that U^ sc \r) provides s.p. state at — 137keV. This 
choice, as we shall see in Sect. IV.A, is essential for a quantitative description of 7 Be(p, 7) 8 B 
radiative capture cross-section. 

The self-consistent average potential U^ sc '(r) (the solid line in the l.h.s. plot) exhibits 
for small r a clear maximum which is absent in the initial potential U(r) . One should also 
notice, that self-consistent potentials U ( - sc \2 + ) and U^ sc \l + ) are different (compare solid 
curves on l.h.s. and r.h.s. of Fig. 1) , in spite of the fact that the initial potential U(r) in 
these states is the same. This clearly shows strong state dependence of calculated average 
fields and two-body effective interactions. For example, whereas the s.p. orbit in the 
g.s. (J 71- = 2±) is at — 137keV , the self-consistent procedure yields this orbit at — 2.538 MeV 
in the first excited state (J n = if) . 

The dotted lines in Fig. 1 show the equivalent s.p. average potentials U^ eq '(r) in J w = 2 + 
and 1 + states for proton lps/2 orbit. For the same a, i?o and Vso parameters as in the initial 
potential U(r), parameter Vq in U^ eq \r) is adjusted in order to reproduce the energy of lps/2 
s.p. orbit in U^ sc \r). Clearly, U^ eq \r) and U^ sc \r) differ strongly in the potential interior. 
However, the potential surface of U^ eq \r) and U^ sc '(r) for both considered J n is very similar. 
This is particularly apparent for J n = 2 + . The root mean squared (r.m.s.) radius of the 
lps/2 orbit is < r 2 > 1/2 = 4.228 fm in £/( sc )(2 + ) and < r 2 > 1/2 = 4.239 fm in U {eq) {2 + ). For 
J w = 1 + , the effective surface diffuseness of U^ sc \r) is even slightly decreasing. Also the 
radial wave functions of the proton lps/2 orbit in U^ sc '(r) and U^ eq \r) are almost identical 
with only a small shift towards the potential interior of the maximum of radial wave function 
in U {sc \r). 

We have no clear indication concerning the position of proton lpi/2 s.p. orbit. Using the 
same U(r) as used to determine U^ sc \r) for lpz/2 s.p. state, we get the lpi/2 proton s.p. orbit 
in U( sc \r) at e Pl/2 = +0.731 MeV in J 71 = 2+ states and at e Pl/2 = +0.311 MeV in J w = 1+ 
states . Consequently, the energy splitting of p^/2 and pi/2 orbitals is state dependent. It 
is Ae = 3.369 MeV, in the initial potential U(r). In U {sc) {r), it equals Ae = 0.868 MeV 
for J n = 2+ and Ae = 2.849 MeV for J w = 1+ states. In J 71 = 3+ states, coupling to the 
continuum does not introduce any correction to the position of lpi/2 s.p. state. Therefore, 
to obtain the energy and the wave function of lpi/2 state we take the equivalent potential 
U^ q \?> + ) to U (sc \3 + ) for the lp 3/2 . 

In general, the surface region of average potential shows weak sensitivity to the self- 
consistent correction. Deviations from this general rule can be seen for weakly-bound many- 
body states having an important admixture of I = and I = 1 neutron s.p. states, i.e., for 
halo configurations in the neutron rich nuclei ||32|| . In this case, the self-consistent correction 
may change significantly the surface of the average potential [|33| . In the above discussion we 
have considered effect of coupling between P and Q subspaces on the average s.p. potential 
for the proton orbits. To exclude the effect of Coulomb barrier, in Fig. 2 we show the initial 
potential U(r), the self-consistent potential U^ sc \r), and the equivalent potential U^(r) 
for neutron lp^/2 s.p. orbit in 8 Li where, for the sake of argument, 1^3/2 neutron s.p. orbit in 
(r) is at —20 keV. Parameters of the initial potential in this case are: Vq = —33.4 MeV, 
Ro = 2.4 fm, a = 0.52 fm and Vso — —4 MeV. Parameters of the residual interaction @ 
are the same as used in Fig. 1. One may notice that the self-consistent correction to the 
average potential is now much stronger than in 8 B and, in particular, it changes strongly 
the surface of [7^ c ^(r), as expected for weakly-bound halo nuclei 
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Fig. 3 shows the dependence of U^ sc '(r) on the parameter (1 — a) of the spin-exchange 
component of the two-body residual interaction ([|). The self-consistent potentials for differ- 
ent values of (1— a) have been determined for the lps/2 s.p. proton orbital (e P3/2 = — 137 keV) 
in J 71 = 2 + . Parameters of corresponding initial potentials U(r) are shown in Table |I[ We 
see that U^ sc '(r) strongly depends on the the spin-exchange term of the residual force @. 
The largest correction to the initial potential U(r) is obtained for small values of (1 — a), i.e., 
when approaching the Wigner force limit. This dependence of U^ sc \r) on the spin-exchange 
component of the residual interaction has a sensible effect on the width of resonances, as we 
shall see below. 



B. Spectrum of 8 B 

In this section and in Sect. III.C, we shall present the detailed calculations of spectra for 
mirror nuclei : 8 B and 8 Li. The SMEC results depend on the following ingredients : (i) the 
nucleon - nucleon interaction in Q subspace, (ii) the residual coupling of Q and P subspaces, 
(hi) the self-consistent average s.p. potential which generates the radial formfactor for s.p. 
bound wave functions and s.p. resonances, and (iv) the cutoff radius for s.p. resonances. 



The freedom of choosing asymptotic conditions in solving eqs. (10) means that a zero of 
excitation energy scale can be fixed arbitrarily. In the following, we choose the zero on this 
energy scale by requiring, independently for any SMEC calculation shown in Fig. 4 and in 
Tables [Q| and |T|, that the lowest experimentally known resonance (J n = if state for 8 B) 
with respect to the proton emission threshold has its energy equal to its experimental value. 
The same convention is used in the calculation of the 7 Be(p, 7) 8 B capture cross section in 
Sect. IV.A. 

Fig. 4 compares SM energy spectrum for T = 1 states of 8 B (on the l.h.s.) calculated 
in the p - shell using the CK interaction [|3T| , with those obtained in the SMEC in different 



approximations. The experimental data are plotted on the r.h.s. of this figure. For the spin- 
exchange parameter, we take (I — a) = 0.27, which is a standard value resulting from a fit to 
the giant dipole resonance in 16 O p4]J34]j35|1 . Different sets of coupling matrix elements are 
taken into account. In the column denoted 'SMEC, we include coupling matrix elements 
between the J w = 3/2~ g.s. wave function of 7 Be and all considered states in 8 B. 'SMECl' 
and 'SMECl*' columns show the results when also coupling matrix elements between the 
J 77 = l/2~ first excited state in 7 Be and states of 8 B are included. In 'SMECl', this 
J 71 = l/2~ state is at the energy (E* = 1.07 MeV) predicted by SM, which differs from the 
experimental energy (E* = 0.429 MeV) of this level. For that reason, in the column denoted 
'SMECl*', we show results of calculations where energy of J 71 = l/2~ state is put at its 
experimental position, without changing neither the coupling matrix elements of residual 
force nor the effective interaction in Q subspace. 

The iterative procedure to correct U{r) and to include the diagonal part contribution 
of residual interaction has been described in the previous section. The self-consistently 
determined s.p. potential U^ sc '(r) is then used to calculate the radial formfactors of coupling 
matrix elements and the s.p. wave functions. One and the same potential U(r) is used for 
the calculation of self-consistent potentials for all many-body states in 8 B, and for both 
1^3/2 and lpi/2 proton s.p. states. These self-consistent potentials for lps/2 proton s.p. 
state in J 71 = 2 + and 1 + states of 8 B can be seen in Fig. 1. For neutrons, there is no 
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correction from the residual interaction, and the average s.p. potential is chosen such that 
it yields lp?,/2 and lpi/2 neutron orbits at — 13.02 MeV and — 11.16 MeV respectively. These 
values have been deduced from experimental Q-values and SM spectroscopic factors. The 
parameters of initial potential for neutrons are : Vo = —67.595 MeV (for lps/2 state) and 
Vo = —70.942 MeV (for lpi/2 state). The remaining parameters are the same as for protons. 

The r.m.s. radius of neutron lp^n and lp\/2 orbitals is < r 2 >^ 2 ~ 2.5 fm. We have checked 
that this choice of potential parameters for neutrons is not crucial and, e.g., shift of neutron 
s.p. energies by ~ 2 MeV does not influence final results of SMEC for the spectroscopy of 8 B. 
Supplementary informations concerning results shown in Fig. 4 can be found in Table |T|. 

The dependence of the 8 B energy spectrum on the strength VyP of the residual interaction 
(0) is shown in Table 0. The SMEC calculations correspond to the standard value of the 
spin-exchange parameter: (1 — a) = 0.27. The calculated energy spectrum depends in a 
complicated way on the strength parameter Vilp. The optimal value of VfjP, in particular in 
what concerns the width of resonance states which provide a most sensible test of the wave 
functions and the values of matrix elements of residual coupling, should be searched in the 
interval between 600 and 700 MeV-fm 3 . In the following, in all SMEC calculations presented 
in this work we use Vy} = 650 MeV-fm 3 , and discuss in some details the dependence of 8 B 
and 8 Li spectra on the spin-exchange term in the residual interaction ([3|). 

The spectrum of 8 B is relatively insensitive to certain approximations in SMEC. Ground 
state energy relative to the proton emission threshold is reasonably well reproduced by the 
SMEC. The position of g.s. with respect to the energy threshold for proton emission changes 
by about 100 keV due to the inclusion of coupling to the g.s. of 7 Be (compare columns 
denoted 'SM' and 'SMEC in Fig. 4 and in Table |IIJ). The coupling to the excited 7 Be* 
nucleus is relatively unimportant (see columns denoted 'SMEC1' and 'SMEC1*' in Fig. 4 and 
in Table [Tip. Calculated width of J 71 " = if depends weakly on chosen couplings (compare 
'SMEC with 'SMEC1' or 'SMEC1*' columns in Table 0) and for the interaction with 
(1— at) = 0.27, one obtains approximately half of the experimental width. Width for J 71 = if 
state depends sensitively on the proportion of direct and spin-exchange terms (see Table fV\j 
in the residual interaction, and the agreement with the experimental value improves when 
approaching the limit of Wigner force (a = 1). In this limit, the residual interaction Vyi in 



(0) is compatible with SU(4) supermultiplet symmetry [SB]. It is well known that this limit 
is badly broken in majority of heavier nuclei mainly due to the increasing importance of spin- 
orbit coupling. However, in p-shell nuclei, SU(4) is still a good approximation |3"7|J3"8(]. The 



Gamow- Teller transitions have been studied using the SM calculations with CK interaction 
showing that the supermultiplet symmetry of SU(4) is well preserved [|38|1 . Also realistic SM 



calculations have confirmed that g.s. wave functions have very small admixture of SU(4) - 
breaking configuration , both at the beginning and at the end of p - shell [31,37,39] . For that 
reason, it is encouraging that the residual two-body interaction between Q and P, which is 
used in SMEC calculations for low-lying many body resonances in 8 B, is consistent with the 
approximate validity of SU(4) scheme found in the SM. 

It is instructive to see the energy dependence of J 71 " = if lowest eigenvalue of the effective 
Hamiltonian flT2Tj in 8 B (see Fig. 5) . The same interaction and potentials as in Fig. 4 are 
used. Results shown with the solid line include the coupling of J 71 = if state to the 
channel wave function obtained by the coupling of J w = 3/2~ g.s. wave function of 7 Be with 
one proton scattering wave function. The dashed line shows results obtained by including 



10 



couplings both to the g.s. and to the J n = l/2~ first excited state of 7 Be at the energy given 
by the SM calculation. Finally, the dotted line shows the results when this first excited state 
is placed at the experimental position (E* = 0.429 Me V). The real part of the eigenvalue 
depends on the total energy of the system. For the center of mass (cm.) energies less 
than about 2 MeV, Er is a decreasing function of energy and for higher energies it becomes 
an increasing function of energy. The minimum of Er(E) depends slightly on whether the 
coupling to the J n = 1/2" first excited state of 7 Be is included. On the contrary, the 
imaginary part of the eigenvalue increases monotonously with energy from the threshold 
and saturates at higher energies. Different approximations concerning the couplings have 
an influence on Er and are practically insignificant for Tr. 

SM energy of the first J n = 3f level is too low as compared to the experimental value 
(see Fig. 4 and Table [TV] ) . The coupling to the continuum cannot correct for this deficiency 
of SM. The width of 3^ state differs by at least a factor 5 from the experimental data 
and here, again, the agreement between experiment and calculations improves when a — > 1 
(see Table |TVj ). There may be several reason for the observed discrepancy between SMEC 
calculations and the experimental data for this state. Firstly, as mentioned already, SM is 
not well describing energy of this state. (The CK interaction is somewhat less successful in 
describing spectra of light p-shell nuclei.) Secondly, the wave function of experimental 3f 
state is most probably overlapping with the cluster configuration [ 3 He — 5 Li], which cannot be 
adequately described in p-space SM calculations. It is known in mirror nucleus 8 Li, that the 
3~± state is strongly influenced by the mirror cluster configuration [ 3 H — 5 He] [^TJ . Moreover, 



energy of experimental 3^ state lies above the threshold for three-particle decay: 8 B — > 
[ 3 He — p — 4 He]. This decay channel which surely contributes to the experimental width, 
cannot be taken into account in the one-particle approximation for scattering continuum. 
Also correlations which could be generated in Q subspace by the coupling to this three- 
particle decay channel are outside of considered model space. Thirdly, 3f state in 8 B cannot 
couple to the first excited state 1/2 - of 7 Be but could couple to higher excited, particle- 
unstable states such as 7/2~. These couplings cannot be treated fully consistently in SMEC, 
because the closeness condition (jl]) implies that states in (N — 1) - nucleus must be stable 
with respect to the particle emission. 



C. Spectrum of 8 Li 

Table compares SM energy spectrum with those obtained in the SMEC for T = 1 states 



of Li. SM calculations in the p - shell (Q subspace) are done using CK interaction |3T 
like in the case of 8 B. Parameters of the residual interaction (Rf) are also the same as given 



in caption of Table III . Parameters of U (r) are given in Table |. U(r) is such that the 



neutron lps/2 s.p. state in U^ sc \r) is bound by 2.033 MeV in J 71 " = 2f state, similarly as the 
experimental g.s. of 8 Li. The lpi/2 neutron state is then bound by 1.109 MeV. The r.m.s. 

of IP3/2 neutron s.p. orbit in (r) is < r 2 > 1/2 ~ 3.474 fm. The same initial potential is 
used for calculation of self-consistent average potentials for all other many-body states in 
8 Li and for both lp^/2 and lpi/2 neutron s.p. states. 

Experimentally, 3i~ state in 8 Li lies above the neutron emission threshold and we fix zero 
on the energy scale by requiring that for all different examples of SMEC calculations shown 
in Tables |7| and [V| the energy position of this lowest resonance with respect to the neutron 
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emission threshold corresponds to the experimental value. The same convention is used also 
in the calculation of the 7 Li(n, 7) 8 Li capture cross section in Sect. IV. B. 

Coupling to the continuum induces the renormalization of spin-orbit interaction. The 
energy splitting of 1^3/2 and lpi/2 orbitals which equals : Ae = 2.302 MeV in U(r), becomes 
Ae = 0.924 MeV in U {sc) (2 + ), and Ae = 3.078 MeV in U {sc) (l + ) . One should notice that 
these splittings are different from analogous splittings in 8 B (see the discussion in Sect 3.B). 
Hence, the coupling to the continuum breaks explicitly the mirror symmetry of spectra and 
wave functions of SM. 

The dependence of self-consistent correction to the average potential on isospin of s.p. 
states, induces a salient dependence of neutron s.p. potential on proton s.p. potential. This 
dependence is very weak in the studied cases of 8 Li and 8 B. For protons, U(r) is such that 
it yields 1^3/2 orbit at —14.8 MeV and lpi/2 orbit at —13.9 MeV. These values have been 
deduced, similarly as for 8 B, from the experimental Q- values and SM spectroscopic factors. 
The parameters of U(r) are V = —73.066 MeV (for lp^/2 state) and V = —78.058 MeV (for 
lpi/2 state), and remaining parameters are the same as for neutrons. 

In the column denoted 'SMEC in Table [V], we include the coupling between J n = 3/2~ 
g.s. of 7 Li and all considered states of 8 Li. 'SMECT and 'SMECl*' columns show the results 
when also the coupling between J 71 = 1/2" first excited state in 7 Li and states of 8 Li are 
included. In 'SMECl', this J w = 1/2" state is at the SM energy : E* = 1.07 MeV. This 
energy value differs from the experimental value : E* = 0.467 MeV. For that reason, in 
'SMECl*' column we show results where the energy of J n = l/2~ state is placed at its 
experimental position. The coupling matrix elements of the residual force (||) and the CK 
interaction in Q subspace remain unchanged in 'SMECl' and 'SMECl*' calculations. 

The position of g.s. with respect to the energy threshold for neutron emission is weakly 
modified by the inclusion of coupling to the continuum. The energy shift in this case is 
smaller than the shift caused by the coupling to J 71 " = l/2~ excited state of 7 Li*. In general, 
the coupling of 8 Li to 7 Li* plays more important role than the analogous coupling of 8 B 
to 7 Be* . Ground state energy with respect to the neutron emission threshold is not well 
reproduced by SMEC calculations. Energies of lowest three states : J n = 2f, if, 3f in 
SM calculation, are too much compressed as compared to the data. On the other hand, the 
energy splitting between 3+ state and the second excited state if is well reproduced. The 
width of 3^ state does not depend on the chosen coupling matrix elements and is by factor 
~ 10 too small as compared to the data. The width of J 71 " = if state depends also weakly 
on the coupling matrix elements (see columns denoted by 'SMECl' and 'SMECl*') and is 
comparable to the experimental decay width. 

Table [VI| compares SM energy spectrum for T = 1 states of 8 Li (on the l.h.s.) with those 
obtained in the SMEC for different parameters (1 — a) of the residual interaction (|3]). A 
satisfactory agreement with the data for the width of if state is obtained when approaching 
the limit of Wigner force (a = 1). Also the width of 3f state improves largely and differs 
only by factor 3 from the data. We consider that this agreement is satisfactory in view of the 
obvious limitations : (i) SM is not well describing energy of this state, (ii) the wave function 
of the 3f state is strongly overlapping with the cluster configuration [ 3 H — 5 He] which cannot 
be reliably described in p - shell calculations |4"0" |, and (iii) 3f state cannot couple to the 



first excited state 1/2 of 7 Li but could couple to higher excited states, such as 7/2 , which 
are particle unstable and cannot be included in the SMEC model. One should notice that 
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SMEC describes better 3J in 8 Li than in 8 B due to the fact that the three-particle decay 
channel in 8 Li is closed. 

Contrary to the J n = 3^ state, the J 71 = state in 8 Li is well reproduced, in particular, 
when the parameter of spin-exchange component in the residual interaction (D) approaches 
the limit of the Wigner force. It is interesting to see the energy dependence of second 
eigenvalue of the effective Hamiltonian (0). The results shown with the solid line in Fig. 
6 include the coupling of J n = in 8 Li to the J n = 3/2~ g.s. of 7 Li. One neutron is in the 
continuum. The dashed line presents results obtained by including couplings both to the g.s. 
and to the J 71 = l/2~ first excited state of 7 Li at the energy predicted by the SM. Finally, 
the dotted line shows the results when this first excited state is placed at the experimental 
position E* = 0.467 MeV. Contrary to the the case of J 77 = if lowest eigenvalue in 8 B, now 
also the imaginary part Tr of the eigenvalue depends on the chosen couplings. 



D. The quadrupole moment of 8 B and 8 Li 

The quadrupole moment < Q > of 8 B provides a useful test of the SMEC wave function, 
in particular of its radial part. It can be calculated following the approach of Carchidi et 
al ||42|| , which has been recently applied by Brown at al [T5[] to the case of 8 B and 8 Li. The 



quadrupole moment (~ < r 2 Y^ >), is obtained by summing over products of many-body 
matrix elements and s.p. matrix elements [ [42]] : 

< J^Y^U >= J2 4 < M\[at t a f J^U ><j,t z \\r 2 Y^\\j',t z > , (16) 

where the sum over t z runs over protons and neutrons. Inserting the complete set of states of 
the A = 7 system between the operators af tz and ay tz , one can write fllBD as a sum over all p 
- shell states of A = 7 system [[43] with t z = p,n , for 7 Be and 7 B respectively. Consequently, 
one can express the matrix element fllBD in terms of corresponding spectroscopic amplitudes, 
effective charges, geometrical factors and radial matrix elements [P~5| . Different radial matrix 



elements which enter the final expression for < Q > (see eq. (20) of Ref. |15 ) correspond 
to the terms coming from the g.s. of 7 Be (the so-called valence term : < r 2 >^ p * ) ), and 
the excited states in 7 Be (the p - core proton term : < r 2 >^pl c ) and in 7 B (the p - core 
neutron term : < r 2 >p_ c ), respectively [|TJ,0. For our radial formfactors, calculated 
using the self-consistently determined average potential for the spin-exchange term with 
(1 — a) = 0.05, the radial matrix elements are : < r 2 >^= 17.9 fm 2 , < r 2 >pl c — 7.65 fm 2 
and < r 2 >pl c — 6.03 fm 2 . The position of p - core proton and neutron states have been 
deduced from the experimental Q - values. With the above values of radial matrix elements 
, one finds : < Q >= 6.99 e fm 2 , in the good agreement with the experimental value E4 



: < Q >= 6.83 ± 0.21 e fm 2 . This theoretical value, which is dominated by the valence 
term, has been obtained assuming the effective charges : e p = 1.35 , e n = 0.35, and the 
SM spectroscopic factors for the CK interaction. The analogous calculation in 8 Li yields : 
< Q >= 2.78 e fm 2 , close to the experimental value [[44J : < Q >= 3.27 ±0.06 e fm 2 . In this 
case, the calculated radial matrix elements are : < r 2 >^= 12 A fm 2 , < r 2 >p1^ c = 7.45 fm 2 

and < r 2 >^l c = 5.9 fm 2 , respectively. Even though the agreement between those theoretical 
estimates of < Q > in 8 B and 8 Li and the corresponding experimental values is encouraging, 
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nevertheless it should not be overstated in view of some uncertainties concerning values of 
the effective charges, the spectroscopic factors and the radial matrix elements for the core 



terms < r 2 >'fl c and < r 2 >i_, . 



„2 ^(n) 



IV. RADIATIVE CAPTURE PROCESSES INVOLVING 8 B AND 8 LI NUCLEI IN 

THE FINAL STATE 

A. The 7 Be(p, 7) 8 B reaction 

The f3 + decay of 8 B , which is formed by the reaction 7 Be(p, 7) 8 B at the cm. energy of 
about 20 keV, is the main source of high energy solar neutrinos. In the absence of agreement 
between different experimental data for this reaction and in view of the disagreement between 
different measurements of solar neutrinos, the input of SSM |l[ should be compared with 
the theoretical values. The cross-section for the reaction 7 ~Be(p, 7) 8 B remains the main 
uncertainty in the input of the SSM. 

In the SMEC, the initial wave for the system [ 7 Be + p] is : 



%.(r) 



5> 

laja 



la ^UjJ 



r 



[Y L x xT x X h 



(17) 



and the final wave for the system [ 8 B]jtt =2 + is: 

, J f 



kjb 



[Y tb x xT x X 



m f 



(18) 



I t and s denote the spin of target nucleus and incoming proton, respectively. A^^ f is the 

j 

coefficient of fractional parentage and Ui j is the s.p. wave in the many-particle state Jf . 
With the wave ^(r) and ^ f{r) , we calculate the transition amplitudes: 



T EC = C(EC)i L J f l b j b j a < C5J f m f \ >< l b 0£0 \ l a > 
yiW{ 3b I t CJ l J f3a )W{l b sC 3a3b l a )iy: h0b 

for El and E2 and : 

rpMl _ i l ^jj, N j f < ISJfTTlf | JiUli > 
xlw(j b I t lJi] Jfja)]a 3 b 

^{p t + ^)UaW{l b slj a ;j b l a ) + (-iy b -i a 2fi a ssW(sl b l 3a ; j b s 

+»t{-i) J !- Ji iti t w{i t j h ij^^ 



(19) 



(20) 



for Ml transitions, respectively. In the above formula, 5 = rrii — rrif, a = \J1a + 1, a = 



I a[a 



+ 1) and : 
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t&M = J u ^ r ^iX dr ■ ( 21 ) 
The radiative capture cross section can then be expressed as: 

K1,M1 167T/fc 7 \ 3 / /i\/e 2 \ 1 1 V~~^ I rriEX.Ml 1 2 



a 



9 \k p J \hcj \hcj 2s + 1 2I t 



•■-^(^(eJsttsettEit-p (23, 

where It and s denote the spin of target nucleus and the spin of incoming proton, respec- 
tively, n stands for the reduced mass of the system. 

Fig. 7 shows the calculated multipole contributions to the total capture cross section 
as a function of cm. energy for different parameters of the residual interaction (|3|). In the 
upper part of the figure, the calculation is done for the same parameters of the residual 
interaction as used in the calculations of spectra shown in Fig. 4. In the lower part, the 
calculation is done for (1 — a) = 0.05, close to the Wigner force limit. The spectrum in this 
case can be seen in Table [TV]. Parameters of the initial potential U(r) in these two cases 
can be read from Table |. The zero of the excitation energy scale is chosen as described 
in Sect. III.B. With this convention, the photon energy is given by the difference of cm. 
energy of [p® 7 Be]j+ system and the experimental energy of the 2f g.s. of 8 B. As can 
be realized from Fig. 7, the El and E2 contributions as well as the total cross-section are 
insensitive to the amount of spin-exchange in the residual force. On the contrary, the Ml 
contribution and particularly its resonant part, are strongly dependent on a. Ratio of E2 
and El contributions at the position of if resonance is 8.15-10 -4 or 7.72-10~ 4 depending 
on whether the spin-exchange parameter equals 0.27 (the upper part of Fig. 7) or 0.05 (the 
lower part of Fig. 7). The experimentally deduced value for this ratio 6.7jlf;| • 10~ 4 |45| is 
consistent with our finding but does not allow to distinguish between different spin-exchange 
parameters. 

El component provides the main contribution to the total capture cross-section in the 
reaction 7 Be(p, 7) 8 B . This non-resonant contribution is a good measure of the spatial 
extension of 2\ wave function, which in turn is determined by the extension of the proton 
1^3/2 orbit. It is essential for the calculated cross-section that the lp^/2 proton orbit in the 
self-consistent average potential is bound by 137 keV. Modification of this value by different 
choice of the depth parameter Vq in U(r), introduces the change in S E1 which can be much 
larger larger than the change due to uncertainties in the potential radius Ro or its surface 
diffuseness a. 

Fig. 8 shows the total S - factor as a function of the cm. energy. The SMEC results 
correspond to a = 0.95. Different multipole contributions to the total cross section for 
this parameter of the residual interaction (§) have been shown in the lower part of Fig. 7. 
Together with SMEC results for the 5-factor, we show experimental data |6|]7]]. The low 
energy dependence of S(E) can be fitted by fll3~l : 



S(E) = S(0) exp(aE + (]E 2 ) . (24) 

In our case, the fit of S(E) using ( p4|) in the range of cm. energies up to 100 keV yields 
5(0) = 19.594 eV-b, a = -1.544 MeV -1 , $ = 6.468 Me V~ 2 . As compared to the similar 
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we find smaller S(0) 
Recent experimental 



fit of S - factor calculated in the Generator Coordinate Method 
parameter and slightly bigger values for the parameters at and (3. 
determination of S-factor p] yields similar to ours low value for S(0). 

The ratio of Ml and El contributions for a = 0.95 is: a Ml /a El = 1.43-10" 3 , 2.65-10" 3 
and 1.90- 10~ 2 at 20, 100 and 500 keV, respectively. The resonant part of Ml transitions 



yields S 



Ml 



20.52 eV-b at the 1 1 resonance energy. This Ml-contribution to the astro- 



physical S-factor decreases fast and becomes S M1 = 3.65T0 -1 , 4.74-10" 2 , 2.72-10 -2 eV-b at 
cm. energies 500, 100, 20 keV, respectively. At the position of if resonance, the calculated 
S - factor (S = 40.67 eV-b) is smaller than measured by Filippone et al. [[J . This value , 
which is dominated by the Ml-contribution, is proportional to the square of spectroscopic 
amplitude of p-states, which for the CK interaction is —0.352 and 0.567 for py% and p%/2 
respectively. Similar small values of spectroscopic amplitudes are obtained for Kumar |46 
and PTBME |47H interactions (see also 
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The E2 contribution to the astrophysical factor was recently measured by Kikuchi et 

0.0^g ( ' +3 ' 6 ' ) meV-b in the energy intervals 



al. who finds: S E2 ~ 0.0] 



° 8 meV-b and S E2 



from 1.25 to 1.5 MeV and from 1.5 to 1.75 MeV respectively. SMEC gives for this quantity 
52 — 53 meV-b and 53 — 71 meV-b, in these two energy intervals respectively. Similar values 
for S E1 have been been found by Typel and Baur |49[|. These values are by factor 10 larger 



than those determined by Kikuchi et al. [48] what remains a puzzle. 



B. The 7 Li(n, 7) 8 Li reaction 

This mirror reaction to the above considered capture reaction: 7 Be(p, 7) 8 B , together 
with a simultaneous description of energy spectra and particle decay widths of 8 B and 8 Li, 
provides another stringent test for SMEC calculations. The SM interaction and SM many- 
body wave functions (e.g. the spectroscopic amplitudes) are identical in both cases. The 
self-consistent one-body potentials which take into account residual coupling of Q and P 
subspaces and which determine the radial formfactors of s.p. wave functions used in the 
calculation of matrix elements of the residual interaction (|3]), are optimized in the same way 
in 8 B and in 8 Li. Finally, the parameters of direct and spin-exchange terms in the residual 
interaction (H) are also the same, so the modification of coupling matrix elements in 8 B and 
8 Li is solely due to the different radial shape of s.p. wave functions in the corresponding 
self-consistent potentials for different J n of many-body states. In the case of neutrons, 
the integral in Eq. ( pip is sensitive to the nuclear interior even in the low energy limit. 
From elastic scattering of neutrons the scattering lengths as, where S is the channel spin 
(S = s + I t , c.f. Eqs. (|I7p and flT8|)), are known to be |5T|: a x = 0.87 ± 0.07fm and 



a.2 = —3.63 ± 0.05 fm. So for the s-wave in the initial channel we adopted procedure of 
Barker [|5l[ of readjustment of appropriate s-wave scattering potentials in order to reproduce 
these experimental values of scattering lengths. 

In Fig. 9 we show different multipole contributions to the total capture cross section as a 
function of cm. energy for different parameters a of the residual interaction (|3p. The same 
parameters have been used for the mirror reaction: 7 Be(p, 7) 8 B (see Fig. 7). In the lower 
part of Fig. 9, the calculation is done for the spin-exchange parameter equal (1 — a) = 0.05. 
The corresponding parameters of U(r) can be found in Table |. The choice of zero on the 
excitation energy scale is the same as described in Sect. III.C. With this convention, the 
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photon energy in Fig. 9 is given by the difference of cm. energy of [n® 7 Li]j+ system and 

i 

the experimental energy of the 2f g.s. of 8 Li . As can be seen in Fig. 9, the total cross- 
section and the El contribution in particular, are insensitive to the amount of spin-exchange 
in the residual force. The E2 contribution shows a weak sensitivity to the parameter a in the 
region of 1^ resonance. The Ml contribution and in particular its resonant part, are strongly 
dependent on a. At the thermal neutron energies, Ml contributions for (1 — a) = 0.23 and 
0.05 differ by approximately one order of magnitude. 

Like for the mirror reaction 7 Be(p, 7) 8 B, the dominant contribution to the total capture 
cross-section in 7 Li(n, 7) 8 Li reaction is the El component. Nevertheless, the Ml contribu- 
tion in 7 Li(n, 7) 8 Li is relatively more important, in particular near the 3f resonance. This 
is partially due to the smaller extension of lp^/2 neutron s.p. orbit in the g.s. wave function 
of 8 Li (lps/2 neutron s.p. orbit is bound by 2.033 MeV in the g.s. of 8 Li) as compared to the 
extension of lps/2 proton s.p. orbit in the g.s. of 8 B. This strong binding of lpz/2 neutron 
state in 2 + states of 8 Li has also a direct consequence on calculated radiative capture cross 
section which becomes reduced, mainly its El component. One should also underline that in 
low-energy reaction such as 7 Li(n, 7) 8 Li , neutron penetrates interior region of the potential 
and, therefore, is more sensitive to the ratio of direct and spin-exchange terms in (^) and, 
indirectly, to the modifications of the interior of average potential by the coupling to the 
continuum (see Figs. 1-3). 

Fig. 10 shows the total neutron capture cross-section as a function of the cm. energy. 
The calculation is done for the same parameters of residual interaction (^) as used in Fig. 8 
for the mirror reaction 7 Be(p, 7) 8 B . Together with the SMEC results, we show the data of 



Nagai et al. which measured the cross-section for the 7 - decay to the g.s. of Li. The 
calculation fits well the data at this energy. The low energy dependence of calculated total 
neutron capture cross section can be fitted by |52|: 



v(E) = ^) 1/ \so + s 1 E^ + S2E + ...) , (25) 

where /i n is the reduced mass of neutron and E is the cm. energy in MeV. In the energy 
interval up to 100 keV SMEC results can be well fitted by three-parameter fit: Sq = 11.517, 
Si = —2.145, s 2 = —11.636, when p n is expresed in a.m.u. 

The ratio of Ml and El contribution is a Ml /a El = 1.78-10" 3 , 2.09- 10" 2 and 1.09 at 20, 
100 and 200 keV, respectively. The resonant part of Ml transitions which is overestimated 
in the calculation due to small calculated width for this state, yields the contribution of 
cr M1 = 305.1 ph at the 3^ resonance energy. This contribution decreases fast and becomes 
a M1 = 13.8, 0.426, 9.18-10" 2 fib at cm. energies 200, 100, 20 keV, respectively. 



V. CONCLUSIONS 

In this work we have applied SMEC model for the microscopic description of spectra in 
8 B and 8 Li, and low-energy radiative capture cross sections in mirror reactions: 7 Be(p, 7) 8 B 
and 7 Li(p, 7) 8 Li. The SMEC model, in which realistic SM solutions for (quasi-)bound states 
are coupled to the one-particle scattering continuum, is a development of CSM model |24]J25 



for the description of low excitation energy properties of weakly bound nuclei. For that rea- 
son, we use a realistic SM effective interaction in the Q subspace and introduce a residual 
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zero-range force with the spin-exchange included which couples Q and P subspaces. This 
deliberate choice of interactions implies that the finite-depth potential generating poten- 
tials in P and matrix elements of residual coupling (Q subspace), has to be determined 
self-consistently. The self-consistent iterative procedure yields new state-dependent average 
potentials and consistent with them new renormalized matrix elements of the coupling force. 
These renormalized couplings and average potentials are then consistently used in the calcu- 
lations of spectra and capture cross-sections, i.e., both in Q and P subspaces. What should 
be taken for coupling between bound and scattering states is in principle not known and we 
have decided to use a schematic combination of Wigner and Bartlett forces. Varying the 
parameter of the spin-exchange component for a fixed intensity of the coupling, we came to 
the conclusion that most satisfactory description of experimental data is achieved for small 
contribution of the spin-exchange part, i.e., approaching the limit of pure Wigner force. This 
finding is consistent with the results of SM which strongly suggest an approximate validity 
of SU{4) symmetry in p-shell nuclei [|TiT , |3"T| -|3l]|j . This proves also an intrinsic consistency 
in our model between the effective SM interaction, in our case the CK interaction, and the 
residual coupling between Q and P subspaces. 

Simultaneous study of mirror system and reactions allows for a better understanding of 
the role that play different approximations and parameters in the model. The dependence of 
final results on radius, diffusivity and spin-orbit coupling parameters of the initial potential 
U(r) is not terribly important and they can be taken from any reasonable systematics. 
Coupling to the excited configurations in the (N — 1) - daughter nucleus is somewhat more 
important in 8 Li than in 8 B. However, this coupling depends only on the wave function of 
daughter nucleus in the excited state and is totally insensitive to the exact energy position 
of these excited configurations. On the contrary, the depth of U(r) has to be carefully 
adjusted so that the energy of s.p. orbit(s) in U^ sc >(r) involved in the systems: [n®(iV — 1)] 
and [p(£)(N — 1)], reflects the binding of g. s. in the nucleus N. This is very important for 
any quantitative analysis of the reaction cross-section. In the studied cases of 8 B and 8 Li 
, the correct identification of this s.p. orbit and hence the determination of an appropriate 
depth parameter in U(r) is simple because the SM spectroscopic factor in 2f g.s. is close to 1. 
Different binding of 8 B and 8 Li leads in SMEC model to different U^ sc \r) for corresponding 
many-body states in mirror nuclei. This in turn causes: (i) the breaking of initial mirror 
symmetry in SM spectra of these nuclei and (ii) the different effective spin-orbit splitting in 
self-consistent potentials for corresponding states in mirror nuclei. 

SMEC model in the present form allows to describe microscopically the coupling to 
one-nucleon continuum. More complicated decay channels like, e.g., those involving the a 
particle in the continuum or more than two particles in the asymptotic states are beyond 
the scope of this model. It is encouraging, however, that these possible shortcomings in the 
description of decay channels, as we have shown on the example of 3^ resonances, are so 
unambiguously reflected in the calculated decay width for these states. In general, the decay 
width is particularly sensitive to the details of the SM wave functions involved and to the 
values of matrix elements of residual coupling so they provide a sensible test of the quality 
of SMEC wave functions and/or approximations involved. 

The overall agreement between experimental data and SMEC calculations for studied 
nuclei is good proving the internal consistency of model assumptions and parameters. The 
astrophysical factor S(0) for the reaction 7 Be(p, 7) 8 B equals 19.424 eV-b which is close to 
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the values reported by Filippone et al. @ and Hammache et al. pf but differs from many 



earlier experiments as well as from many theoretical analysis [12,14|,l5| . The low-energy 



dependence of S(E) is slightly different from the one found by Descouvement and Baye |T3 
which is sometimes used in extrapolating the experimental results to the astrophysically 
relevant region. The calculated ratio of S E1 and S E2 contributions in the region of if 
resonance also agrees well with the data of Davids et al. |j45| . Surprisingly, at higher energies 
(E > 1.25 Me V) the calculated E2 contribution differs by factor 10 from the reported value 



of Kikuchi et al. ||48|| . The results of Refs. and |^8| seem to us incompatible with each 
other. 

The present studies have shown that SMEC results depend sensitively on very small 
number of parameters. Some of them, like the parametrization of the residual interaction 
which couples states in Q and P subspaces, has been established in the present work. The 
others, related to the energy of s.p. states which determine the radial wave function of many- 
body states, are bound by the SM spectroscopic factors and experimental binding energy in 
studied nuclei. This gives us a confidence that SMEC can have large predictive power when 
applied to other p - shell nuclei and to other capture cross-sections in this region. 
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APPENDIX A: THE COUPLING POTENTIAL AND THE SOURCE TERM 

To solve the coupled channel equations in P subspace, one has to calculate the matrix 
elements V4/? 7 «5 of the residual interaction (|) between states in Q and P subspaces. For the 
zero-range force including the spin exchanges, which was used in this work, we have: 



V^iCri, r 2) = Z a(3~,6^( r l ~ r 2) 



where 



7 l 



47T 



[a - b)M L l l5 + (1 - 5 TaT0 5 TlTS )(a + b)M^ s 



(Al) 



(A2) 



and angular two-body matrix element with isospin T can be expressed as [}41 



■i(i + (-i) 



X 



-iyp+3 



L 2 

11, 11, 

<j a --jfs-\L0 >< 3l --j s -\IA > 



1-f-l 



.L+T+lj+ls 



11, 11, r / ,n 

■<Ja^Jf3^\Ll >< h -J S -\Ll > [1 + ("1) J 



(A3) 



Any symbol containing the hat, like e.g. j, means: j = y / 2j+T. 
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The channel-channel coupling in (Q), ( |T0D is taken into account through the matrix 
elements of the type: 



Vrr' 



r L -i U 

V aPjs(( J t3a)j\\ (44) "K 5 *) IIGWj) 



(A4) 



The A-body matrix element in the above expression, is rewritten as : 

o 



A' 



ja J/3 L 
3l 36 L 
A K 



4& 



A" 



7 J ' I a ^<5 



A 



IIG&r)') 



(A5) 



and the reduced matrix element in this latter formula is then expressed as : 

Jt ja J 



J 2 - 



Jt Ja J , , v K r . \ K 

$ % J f ' ^ J 'll (4 a *J H J t) ' O'all [a]Ay) \\h) ■ 



(A6) 



The diagonal parts of this operator induce corrections which renormalize the s.p. average 
potential U(r) (see Eq. @ ). 

The source term in the inhomogeneous eqs. flT0| ) takes into account the couplings of the 
type : 



w® = - y £L-V£ hB ({J t ja)j\\ (44) -(^) L 



Y. L ' V a/3 7 s(( J t3 a )j\\ 4 (4 K^)^ ' 



(A7) 
(A8) 



Again, the reduced matrix element in this expression is given as a coupled product of the 
two contributions: 



■i 1 ■ \ ? ''■> ) ■ (UKmuw (4(^)T ii*j> 



(A9) 



This operator modifies both real and imaginary parts of eigenvalues of Hqq in Eq. (12), 
but does not change the s.p. average potential U(r). 
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TABLES 



TABLE I. Parameters of the initial potentials U (r) (|5|) used in the calculations of self-consistent 
potentials U^ sc \r) for various parameters (1 — a) of the spin-exchange term in the residual interac- 
tion (^). U^ sc \r) are constructed for various T = 1, positive parity states in 8 B and 8 Li and for dif- 
ferent single particle states: L»3/2 an d lpi/2- Lor all cases the radius of the potential is Rq = 2.4 fm, 
the diffuseness parameter is a = 0.52 fm, and the spin-orbit orbit parameter is Vso = — 4MeV. 
The strength of the residual interaction (||) is V$ = 650MeV-fm 3 , for all considered cases. 



System 


e Pi/2 [MeV] 


1 — a 


V [MeV] 


[p (2) 7 Be] 


-0.137 


0.45 


-42.140 






0.27 


-40.045 






0.05 


-37.660 


[n (8) 7 Li] 


-2.033 


0.45 


-41.683 






0.27 


-39.555 






0.05 


-36.905 
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TABLE II. Dependence of the spectrum of B on the strength V^ 2 of the residual interaction 
@ and for a fixed value of the parameter of the spin-exchange term : (1 — a) = 0.27. For more 
details, see the description in the text. 



Nucleus 


J 77 


Energy [keV] 


Width [keV] 


Energy (exp.) [keV] 


Width (exp.) [keV] 


y(0) 


[MeV-fm 3 ] 


8 B 


2+ 


-443 


— 


-137.5 ± 1.0 


— 




500 






-441 


— 








600 






-356 


— 








650 






-339 


— 








700 






-166 


— 








800 




1+ 


637 


2.9 


637 ±6 


37 ± 5 




500 






637 


3.5 








600 






637 


16.5 








650 






637 


19 








700 






637 


360 








800 




3+ 


1229 


8.8 


2183 ± 30 


350 ± 40 




500 






1 O0 1 

izzl 


y.z 








OUU 






1294 


13.1 








650 






1311 


12.1 








700 






1470 


16.2 








800 




1+ 


2279 


32 


not known 


not known 




500 






2253 


42 








600 






2153 


240.2 








650 






2128 


262 








700 






2300 


763 








800 
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TABLE III. SM energies and SMEC energies and widths vs. experimental ones of 8 B nucleus. 
All units are in keV. The proton separation energy is adjusted in order to reproduce the energy 
of the lowest resonance state. Different labels denote as follows: 'SMEC - only the ground state 
of 7 Be was taken into account in all couplings, 'SMEC1' - coupling to the first excited state 
in 7 Be was included with E* = 1.07 MeV (SM value), 'SMEC1*' - the same with excited 7 Be 
state at E* = 0.429 MeV (experimental value). Parameters of the residual interaction (||) are: 
V$P = 650MeVfm 3 , a = 0.73. The cut-off radius is R cu t = 5fm except for the pi/ 2 s.p. wave 
function in 1+ many body states, which is in the continuum at about 300 keV above the threshold 
and for which larger cut-off was used R cu t = 10 fm. The numbers in parentheses are the widths of 
3± state if this state would be placed at the experimental energy. For more details, see the Table | 
and the description in the text. 



State 


SM 
energy 


SMEC 
energy width 


SMEC1 
energy width 


SMEC1* 
energy width 


experiment 
energy width 


2+ 


-446 


-356 




-334 




-329 




-137.5±1.0 




1+ 


637 


637 


16.5 


637 


15.3 


637 


15.2 


637±6 


37±5 


3+ 


1246 


1294 


13.1 


1237 


12.5 


1241 


12.6 


2183±30 


350±40 








(25.2) 




(25.8) 




(25.8) 






1+ 


2327 


2153 


240.2 


2081 


272.7 


2080 


309.0 


not known 



TABLE IV. The dependence of 8 B spectra on the relative strengths of direct and spin exchange 
parts of the residual interaction (|3|). Analogously to the entry 'SMEC in Table III only ground state 
of 7 Be was taken into account. The entries in this table are labelled by the value of a parameter 
of the residual force. The numbers in parentheses as in Table III. For more details see Table |, the 
caption of Table III and the discussion in the text. 



State 


SM 


a = 


0.55 


a = 


0.73 


a = 


0.95 


experiment 


J* 


energy 


energy 


width 


energy 


width 


energy 


width 


energy width 


2+ 


-446 


-418 




-356 




-320 




-137.5±1.0 


1+ 


637 


637 


12.6 


637 


16.5 


637 


25.9 


637±6 37±5 


3+ 


1246 


1313 


3.0 


1294 


13.1 


1275 


34.9 


2183±30 350±40 








(5.6) 




(25.2) 




(67.4) 




1+ 


2327 


2299 


115.9 


2153 


240.2 


1899 


398.9 


not known 
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TABLE V. The same as in Table III but for 8 Li nucleus. The neutron separation energy is 
adjusted in order to reproduce the energy of the lowest resonance state. Different labels denote as 
follows: 'SMEC - only the ground state of 7 Li was taken into account in all couplings, 'SMEC1' - 
coupling to the first excited state in 7 Li was included with E* = 1.07MeV (SM value), 'SMECl*' 
- the same with excited 7 Li state at E* = 0.467 MeV (experimental value). The numbers in 
parentheses are the widths of 1^ state if this state would be placed at the experimental energy. 
For other informations see Table ffl, the the caption of Table ID and the discussion in the text. 



State 


SM 
energy 


SMEC 
energy width 


SMECl 
energy width 


SMECl* 
energy width 


experiment 
energy width 


2+ 


-1471.1 


-1437.09 




-1346.4 




-1345.6 




-2033.8±0.3 




1+ 


-388.1 


-418.9 




-349.5 




-351.5 




-1053.0±0.1 




3+ 


221.2 


221.2 


3.4 


221.2 


3.5 


221.2 


3.5 


221.2±3.0 


33±6 


1+ 


1301.6 


1065.6 


357.9 


1050.0 


383.1 


1036.7 


464.2 


1176 


« 1000 








(378.1) 




(422.9) 




(506.8) 







TABLE VI. The dependence of Li spectroscopy on the relative strengths of direct and spin 
exchange parts of the residual interaction is presented. Analogously to the entry 'SMEC in Table |V|, 
only g.s. of the 7 Li was taken into account. The entries in this table are labelled by the value of 
a parameter of the residual force. The numbers in parentheses as in Table [V|. For other informations 
see Table [jj the the caption of Table [v| and the discussion in the text. 



State 


SM 
energy 


a = 0.55 
energy width 


a = 0.73 
energy width 


a = 0.95 
energy width 


experiment 
energy width 


2+ 
1+ 
3+ 
1+ 


-1471.1 
-388.1 
221.2 
1301.6 


-1514.7 
-441.7 
221.2 0.8 
1199.2 171.8 
(170.1) 


-1437.9 
-418.9 
221.2 3.4 
1065.6 357.9 
(378.1) 


-1330.3 
-394.4 
221.2 9.4 
790.0 561.1 
(735.3) 


-2033.8±0.3 
-1053.0±0.1 

221.2±3.0 33±6 
1176 1000 
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FIGURES 




FIG. 1. Finite-depth average s.p. potentials used to generate the radial s.p. wave functions for 
bound states and resonances. 

(i) Plot on the l.h.s. shows the initial potential U(r) ^) (the dashed line), the self-consistent 
potential U^ sc \r) (the solid line), i.e., U(r) which is corrected by the coupling to the continuum of 
scattering states, and the equivalent average potential U^ eg \r) (the dotted line) of the Saxon- Woods 
type which yields the proton orbit at the same energy as in the self-consistent potential. The 

residual interaction (|J) parameters are: = 650MeV-fm 3 , a = 0.73. U(r) is chosen in such a 
way that U^ sc \r) yields the 1^3/2 s.p. state at the energy — 137keV, corresponding to the binding 
energy of the 2~f g.s. in 8 B. The correction to the average potential from the residual interaction 
@ is calculated for the 1^3/2 s.p. orbit in 2 + states of 8 B. For more details, see Table | and the 
description in the text. 

(ii) Plot on the r.h.s. shows the same as on l.h.s. but for the 1 + states of 8 B. The initial potentials 
on the r.h.s. and l.h.s. are identical and U^ eq \r) is reproducing the position of the 1^3/2 s.p. orbit 
of U^ sc \r) for 1+ states in 8 B. 
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FIG. 2. The initial potential U(r) of Saxon- Woods type (see @) (the dashed line), the 
self-consistent potential U^ sc \r) (the solid line) and the equivalent Saxon- Woods average potential 
Jj( eg \r) (the dotted line) which yields the neutron lp%/2 orbit at the same energy as found in 
U^ir). The curves show the absolute value of corresponding potentials in the logarithmic scale 
so the characteristic tail for large r can be well seen. U(r) is chosen in such a way that U( sc \r) for 
neutrons obtained from it by including the diagonal correction term from the residual interaction 
(|3]) yields the IP3/2 neutron s.p. state in J w = 2 + states of 8 Li at the energy — 20keV. For more 
details, see Table | and the description in the text. 
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FIG. 3. The self-consistent average potential U^ sc \r) which is obtained by including the cou- 
pling of states in Q and P subspaces with the residual interaction (|3|) , is plotted for two parameters 
of the residual force a = 0.55 (the solid line) and a = 0.95 (the dashed line). The calculations have 
been done for the s.p. orbit an d in J w = 2 + states of 8 B. The initial potential in each case 

has been chosen in such a way that the corresponding self-consistent potential yields the s.p. 
orbit at the same energy of — 137 keV. For more details, see Table | and the description in the text. 
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SM SMEC SMEC1 SMEC1* exp 



FIG. 4. SM with CK interaction and SMEC in different approximations labelled 'SMEC, 
SMEC1' and 'SMEC1*' vs. experimental T = 1 states of 8 B nucleus. The proton threshold energy 
is adjusted to reproduce position of the if first excited state. The shaded regions represent the 
width of resonance states. For the details of the calculation, see the description in the text and in 
the caption of Table ffl. 



30 




> 



0.85 - 
0.4 : 
0.3 - 



^ 0.2 - 



0.1 







_l I I I I 1_ 



_l I I I I I I I I I I I 1_ 



_l I I I I I I I I I I I 1_ 



1 2 3 4 5 
E CM [MeV] 



FIG. 5. Energy dependence of the eigenvalue (both real and imaginary Tr parts) of the 
effective Hamiltonian ( |T2| ) for the 1^ state in 8 B. The solid line corresponds to the coupling to the 
g.s. of 7 Be only. The dashed line corresponds to the inclusion of coupling to the first excited state 



in 'Be at the energy predicted by SM with CK interaction [31|. The dotted line corresponds to the 
inclusion of coupling to the first excited state in 7 Be which is placed at the experimental energy. 
For more details, see the description in the text. 
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FIG. 6. The same as in Fig. 5 but for the state in 8 Li. For more details, see the caption of 
Fig. 5 and the description in the text. 
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FIG. 7. Multipole contributions to the total capture cross section of 7 Be(p, 7) 8 B as a function 
of the center of mass energy. The SMEC calculations have been done for different values of the 
spin-exchange parameter (1 — a) = 0.27 (the upper part of the figure) and (1 — a) = 0.05. For 
other details, see Table | and the discussion in the text. 
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FIG. 8. The astro-physical /S-factor for the reaction 7 Be(p, 7) 8 B is plotted as a function of the 
center of mass energy. The SMEC calculations have been done for the residual interaction (g) the 
spin-exchange parameter (1 — a) = 0.05. The experimental points have been taken from Filippone 
et al. H and Hammache et al. f7|. 
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FIG. 9. Multipole contributions to the capture cross section to the g.s. (J n = 2^) of 8 Li in 
the reaction 7 Li(n, 7) 8 Li are plotted as a function of cm. energy. The SMEC calculations have 
been done with the spin-exchange parameter (1 — a) = 0.27 (the upper part of the figure) and 
(1 — a) = 0.05. For other details see the discussion in the text. 
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FIG. 10. The cross-section for the reaction Li(n, 7) 8 Li is plotted as a function of cm. energy. 
The SMEC calculations have been done for the spin-exchange parameter (1 — a) = 0.05. The 
experimental point is taken from Nagai et al. [19]. 
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